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Abstract. The slow neutron capture process in massive stars (the weak s-process) produces 
most of the s-only isotopes in the mass region 60 < A < 90. The nuclear reaction rates used in 
simulations of this process have a profound effect on the final s-process yields. We generated 
ID stellar models of a 25M star varying the 12 C + 12 C rate by a factor of 10 and calculated 
full nucleosynthesis using the post-processing code PPN. Increasing or decreasing the rate by a 
factor of 10 affects the convective history and nucleosynthesis, and consequently the final yields. 



1. Introduction 

Elements in the solar system are formed from a variety of nucleosynthesis processes in stars, 
such as the slow and rapid neutron-capture processes (the s-process and r-process respectively). 
The s-process signature of the solar system abundances can be split into three components; the 
weak-component, which involves nuclei with 60 < A < 90, the main-component, which involves 
nuclei having atomic mass between 90 < A < 208, and the strong component, which accounts 
for the production of the solar 208 Pb [I]. The s-process site is attributed to massive stars for 
the weak component, and to AGB stars having initial mass between 1 and 3M at solar-like 
metallicity for the main component and at low metallicity for the strong component [2]. Stellar 
evolution models of massive stars can be used to determine information on the conditions within 
stellar interiors and calculate the s-process yields relevant for the weak-component. Any changes 
to the input physics, such as improved laboratory nuclear reaction rates, can affect the evolution 
of the simulated star and consequently affect the s-process yields. Therefore, different nuclear 
reaction rates can be tested for their astrophysical impact. In this work, we will show how 
variations in the rate of the 12 C + 12 C reaction affect the evolution of a 25M star and the 
consequential s-process yields of the star. The motivation for this work originates from nuclear 
physics experiments and theory (see for instance Spillane et al. 2007 [3J and Gasques et al. 



2007 [1]), which are probing low enough energies to investigate the reaction within the Gamow 
window for 12 C+ 12 C fusion. 



2. Carbon burning in massive stars and the s-process 

Carbon-core burning occurs at T = 0.6 — 0.8 GK. The dominating reactions are 12 C( 12 C,a) 20 Ne 
and 12 C( 12 C,p) 23 Na which, combined with the efficient 23 Na(p,a) 20 Ne reaction, leaves 16 and 
20 Ne as the dominant isotopes in the core (the 16 O being a remnant from the previous helium 
burning stage). In addition to core burning, there are typically multiple stages of shell burning 
later in the evolution of the star, where the number of shells differs depending on the initial mass 
of the star. The main neutron source for the s-process in massive stars is 22 Ne by the reaction 
22 Ne(a,n) 25 Mg. 22 Ne is produced in helium burning from the 14 N(a, 7) 18 F(/3 + ) 18 0(o;, 7) 22 Ne 
reaction chain during the supergiant phase. 22 Ne is partially burnt at the end of helium burning 
with a neutron density, n n ~ 10 6 cm -3 (see, for example, Raiteri et al. 1991 [5j). The 22 Ne 
left in the CO-core is then burnt during carbon-shell burning with higher neutron densities 
(n n ~ 10 11 — 10 12 cm -3 ), with the 12 C( 12 C,a) 20 Ne reaction providing the a particles [6]. The 
carbon-shell burning is sensitive to the profile of 12 C after core-carbon burning, which is in turn 
sensitive to the still uncertain 12 C(a,7) 16 rate and the choice of convection physics in the 
models [6j[7j[8j. In addition, the yields are sensitive to changes in reactions involving 22 Ne or its 
formation [6] jS]. Most of the s-process-rich material ejected by the supernova event of a 25M 
star is formed by carbon-shell burning, affecting the weak s-component [T0][TT]. 

Stellar models of a 25M star with metallicity Z = 0.01 were generated using the Geneva 
Stellar Evolution Code [12J, with the nuclear network post-processed using the NuGrid PPN 
tool [13] . The post-processing was computed using the KHAOS cluster at Keele University. 
Variations in the carbon burning rate were chosen in relation to the 'standard' Caughlan and 
Fowler rate (CF88 from now on). Models were generated with the CF88 rate (referred to as 
model C12s in this work), the CF88 rate multiplied by 10 (C12tl0) and divided by 10 (C12dl0). 
The ratio of the p- and a-channels was chosen to be 35%:65% [14J. The additional uncertainty 
associated with this choice of ratio will be investigated in a forth-coming paper. 
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Figure 1. Structure evolution diagrams of 25M models at solar metallicity using different 
12 C burning rates. Here, (a) is model C12tl0, (b) is model C12s and (c) is model C12dl0. The 
shaded regions correspond to convection zones present in the star, with the type of burning 
indicated by labels. 



Differences in the convective history can be seen in each of the different models. Fig. [T] 
shows structure evolution diagrams for the models C12tl0, C12s and C12dl0. In each case the 
convection zones for carbon shell burning at log 10 (time in years until core collapse) ~ 1 to -2 
differ in size and duration. For model C12tl0, the duration of the second carbon-shell burning 
process is much longer and the size of the convection zones is larger. A large fraction of overlap 
between the first shell and the second shell is observed and isotopes produced in the first shell 
by the s-process will be mixed into the second. The presence of overlapping carbon shells was 
previously noted by Chieffl et. al. [15J. For model C12dl0, the second carbon-shell burning 
episode occurs later in time than for the other two models, which is due to the star contracting 
further before it can reach a temperature high enough so that carbon burning is activated in 
the shell. 

Fig. [2] shows the relative abundances of stable isotopes in the second carbon shell of the 
C12dl0 model with respect to the C12s model. Fig. [3] shows the equivalent plot for the C12tl0 
model with respect to the CI 2s model. In fig. [2j a clear signature of a higher neutron density 
is shown, e.g. higher production of r-only species 70 Zn, 76 Ge and lower production of 80 Sr due 
to the 79 Se branching. The reason is that the star in model C12dl0 contracts further, resulting 
in carbon shells that burn at a higher temperature than the C12s model. Therefore the 22 Ne 
neutron source is burned at a higher temperature, increasing the neutron density. The s-process 
efficiency is quite similar in model C12s and C12dl0 (see, for example, the similar production 
of neutron magic isotope 88 Sr). Concerning model C12tl0, in fig. [3J a general increase in the s- 
process efficiency for isotopes with 60 < A < 90 is shown, compared to model C12s. The overlap 
between the first and second convective shells causes the initial distribution of isotopes at the 
start of the second carbon-shell burning to be affected by the products of the first shell. In the 
first shell, neutrons are also efficiently produced by an additional neutron source, 17 0(a,n) 21 Ne, 
where 17 is mostly produced by 16 0(p,7) 17 F(/3 + ) 17 0. 




Figure 2. The relative abundances of stable isotopes up to Niobium between models C12dl0 
and C12s. S-only isotopes are indicated with a red circle. 



In this paper, it has been shown that changes to the carbon burning rate by a factor of 10 in 
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Figure 3. The relative abundances of stable isotopes up to Niobium between models C12tl0 
and C12s. S-only isotopes are indicated with a red circle. 



stellar models of a 25M star significantly affect the s-process yields. The overlap of convective 
carbon shells active at different burning temperatures also has important implications for the 
s-process yields. Further analysis will be conducted in a forth-coming paper. 
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